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Abstract 

We present a chirally symmetric hadronic model for the vector and axial-vector 
correlators in vacuum. The dominant contributions to these correlators come from 
intermediate pions, p and a± mesons which are calculated in one-loop approximation. 
The resulting spectral functions are compared with the data obtained by the ALEPH 
collaboration from r decay. In the vector channel we find good agreement up to q 2 = 
1 GeV 2 , in the transverse axial- vector channel up to q 2 = 2 GeV 2 , corresponding to 
the regimes dominated by two-pion or three-pion decay channels, respectively. The 
longitudinal axial correlator is in almost perfect agreement with the PCAC result. 



1 Introduction 



The investigation of matter under extreme conditions and the modifications of hadron 
properties with density or temperature is one of the main topics in intermediate and 
high-energy nuclear physics. Experimentally, the cleanest information is obtained from 
electromagnetic probes which penetrate the medium almost undisturbed. In the vector 
dominance model of Sakurai |1| electromagnetic processes are mediated by neutral vector 
mesons (p°, to and <ft). Therefore possible medium modifications of vector mesons attract 
particular interest. For example, the dilepton data measured by the CERES collaboration 
in ultra-relativistic nucleus-nucleus collisions seem to give evidence for a change of the 
p-meson properties in hot or dense nuclear matter. The strong enhancement of dilepton 
production at low invariant masses (~ 0.3 — 0.6 GeV), was interpreted by Li et al. || 
as a signature for a dropping p-meson mass, conjectured by Brown and Rho M as a 
consequence of scale invariance and partial restoration of chiral symmetry. On the other 
hand the CERES data could also be explained in an alternative way by a strong broadening 
of the p-meson in the medium due to resonance formation [||, || and medium modifications 
of the pion cloud [7|, || ||. (For a review see 

Despite of the phenomenological success of this viewpoint, the role of chiral symme- 
try is not obvious. An important consequence of chiral symmetry restoration is that the 
vector and axial-vector (isovector) correlators, which differ in vacuum due to spontaneous 
chiral symmetry breaking, have to become identical in the restored phase, i.e. at high 
temperature or density. As a precursor, already at much lower temperatures the correla- 
tors start to mix, i.e. the in-medium vector correlator receives contributions proportional 



to the vacuum axial-vector correlator and vice versa |TT . 

Roughly speaking, the vector correlator is dominated by the p meson, while the axial- 
vector correlator gets important contributions from the pion and from the a\ meson. In 
this sense, in the above-mentioned phenomenological models the medium modification of 
the p-meson could be interpreted as some kind of mixing phenomenon, caused by thermal 
or nuclear pions |L_^, [Tj|. Certainly, it would be much more satisfactory, to start from 
a model which is manifestly chirally symmetric and to study both, the vector correlator 
and the axial-vector correlator simultaneously. Among other things this implies that we 
have to include the ai-meson and to treat it on an equal footing with the p. Of course, to 
obtain a realistic description, it is not sufficient to stay at tree-level, where the p and a± 
are stable particles without width, but we should also include the most important decay 
channels p — > tttt and a\ — > pn. This requires at least a one-loop calculation. The aim of 
the present paper is to construct such a model in vacuum, while the application to finite 
temperatures will be the topic of a later publication. 

In the 1960s the gauged linear a model emerged as a combination of the Gell-Mann- 
Levy a model, which exhibits the essential features of the spontaneously broken global 
SU{2) L x SU{2) R chiral symmetry |TJ]], and Sakurai's idea of vector meson dominance 
0, in which a Yang-Mills gauge theory is constructed with the SU{2) V isospin group as 
a gauge group and the p meson as the gauge boson. Analogously, Lee and Nieh [IIS] have 
built a Lagrangian with local SU(2)l x SU(2)r symmetry and p and a\ mesons as gauge 
bosons. The local (gauge) symmetry is explicitly broken only by the vector meson mass 
terms, leading to the current-field identities, and global chiral symmetry is spontaneously 
broken as in the ordinary a model, leading to the p — a\ mass splitting. In order to render 
the theory renormalizable, 't Hooft proposed to generate the vector meson mass by using 
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the Higgs mechanism rather than by explicitly breaking the local symmetry [ |T6|j . 

However, from a modern point of view, there seems to be no reason why in an effective 
hadronic model chiral symmetry should be a local symmetry, while it is only a global 
symmetry in QCD, the underlying fundamental theory. In fact, the main advantage of 
this concept is that it allows only for a small number of free parameters since all couplings 
of gauge fields are determined by only one coupling constant g [[17]]. Unfortunately, the 
"pure" gauged linear (or non-linear) a model as described above (including the vector- 
meson mass term) gives the wrong phenomenology of p and a x mesons. In the literature 
one finds two strategies to cure this problem. The first one consists of adding terms of 
higher dimension (i.e. with more fields or derivatives) to the Lagrangian, which are still 
invariant under local transformations. One term of this kind was already considered in 
Ref. [[18]. In Ref. [IT] it was shown that this is not sufficient, and therefore one more 
term was added. In a more recent work [IS] this strategy was combined with the second 
one, which consists of introducing further gauge-invariance breaking terms in addition to 
the mass term. Of course, these terms destroy the current field identity, but we do not 
believe that this is a severe objection against this strategy. One could even argue that 
the current-field identity cannot be valid in nature, since it is in contradiction to the large 
four and six pion continuum observed in the vector correlator above the p resonance. In 
this article we will completely abandon the concept of local symmetry. We will show that 
in this way we can reproduce rather the experimentally measured spectral functions of 
the vector and the axial- vector correlators without adding terms of higher dimension. 

The paper is organized as follows. In Sect. ^| we introduce the Lagrangian of our model. 
We begin with a brief review of the properties of the gauged linear a model in Sect. 



Then, in Sect. |2.2| , we relax the constraints imposed by the gauge symmetry, retaining only 
a global chiral symmetry. Our model is complete after including electroweak interactions 
in Sect. j2J3| . In Sect. |3| we construct the propagators of the mesons (ir, a, p and ai) 
in a one-loop approximation. Here special care is taken, not to destroy the symmetries 
of the model by an unsuited regularization scheme. Since we give up the current-field 
identity, the vector and axial-vector currents are not just proportional to the fields, but 
involve additional terms. The corresponding correlation functions are derived in Sect. 
|j. In Sect. |5| we present numerical results for the vector and axial- vector correlator, and 
related observables, such as the electromagnetic form factor of the pion and the weak r 
decay. Our results are compared with experimental data. Finally, in Sect. || we draw 
conclusions. 



2 Linear a models with vector mesons 

In this section we want to construct the chirally symmetric Lagrangian which will be the 
starting point for our analysis. In view of the phenomenological success of the vector 



dominance model [|I|] in the vector channel, the gauged linear a model which is its 
generalization to SU(2)l x SU(2)r chiral symmetry, seems to be a natural candidate. In 
the literature, this model served e.g. as a toy model for the calculation of p and a± meson 
properties in a hot meson gas [pO] . However, as we will discuss, the gauged linear a model 
does not correctly reproduce the phenomenology of p and a\ mesons in vacuum, and we 
will have to modify it. It is nevertheless useful to begin with a brief review of the basics 
of this model because most of the definitions and formulas given here will remain valid 
after our modifications. 
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2.1 The gauged linear <r model 

Defining the four-component field 



•-(?). m 

we can write the Lagrangian of the mesonic part of the linear a model as 

C„ = l -d^ ■ - ^ $ • $ - ^ ($ ■ $) 2 . (2) 

For /i 2 < chiral symmetry is spontaneously broken, yielding massless pions. To obtain 
massive pions, we add a small explicit symmetry breaking term: 

£ S B = CCr . (3) 

Let T and T 5 be 4 x 4 matrices acting in the (er, 7?) space, with Tj generating ordinary 
isospin rotation around the 7Tj axis and Tf inducing rotations in the o — 7r< plane. They 
fulfill the commutation relations 

[Tj , 3>] = ze iife T fc , [Tj , Tj 5 ] = ze iife T fc 5 , [Tf , Tf ] = ie ijk T k , (4) 

and are normalized to 

trTT j = 2^, trTT/ = 0, tr T*T? = 2<% . (5) 

One can easily see that these operators generate an SU (2) l x SU (2)r group by inspecting 
the commutation relations of the combinations 

7f = ±(T i -Tf), T l R = ~(T i + Tf). (6) 
The Lagrangian (0) is invariant under global SU(2)l x SU(2)r transformations U, 

$ _> u$ , U = e iS - f+i P- f5 . (7) 
To make it invariant also under local transformations U(x), we introduce gauge fields , 

Y fl = p fl .f + a lfl -f 5 , (8) 
and replace the ordinary derivatives in Eq. (0) by "covariant derivatives", 

T>^ behaves under local transformations in the same way as the field $, provided the 
gauge fields transform as 

Y^UY^+'-Ud^UK (10) 
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The field strength tensor of the gauge fields is defined in the usual way as 

V = d^Y v - d v Y„ + ig [F M , Y v ] =: p^ ■ f + a lfiu ■ f 5 , (11) 

with 

/V = <9 M pl, ~~ ®vPn + 9 P/i x Pu + ga lfl x a lu , 
ai^u = d^aiv - d v a Xil + ga lfl x p u + g p M x ai„ , (12) 

and transforms as Y^ v — > UY^W . Therefore the Yang-Mills term 

£y = -itr VF"" = -i (/v • p"" + a lM „ • aD (13) 

is obviously gauge invariant. 

Finally we add a mass term for the vector mesons, which is invariant under global 
transformations, but breaks gauge invariance: 

2 2 
Ttl Ttl 

£ mo = ?f tr Y,Y» = -f (p, • p" + a M • af ) (14) 
Combining everything, we obtain the final Lagrangian: 

^gauged XV* + £sb + £y + £ mo (15) 

As a next step we will briefly describe how this Lagrangian is usually treated. Due 
to spontaneous symmetry breaking the a field has a non-vanishing expectation value, so 
that it is convenient to redefine the a field, 

Oncw = Oold - cr , $ ncw = $ old - fejj = $ i d - $ , (16) 

and to adjust <7o such that the expectation value of the new a field vanishes (i.e. a is 
the expectation value of the old a field) . At tree level, this is achieved by minimizing the 
potential, i.e. by solving the equation 

A 2 erg + p 2 a - c = . (17) 

After the shift of the a field the masses of pion and a and the masses of p and a\ are no 
longer degenerate: 

2 2 i \ 2 2 2 2 , o\2 2 2 2 2 2,2 2 /i o\ 

"V = A* + * 0b » m CT = p+3Aa , m p = m , m ai = m + # ct . (18) 

Eq. (|T7| ) can now be written as m^ o"o — c = 0, i.e. for c = the pions are massless, as 
required by the Goldstone theorem. 

In addition to the mass splitting, the shift of the a field generates many new vertices, 
in particular a ir — ai vertex generated by the interaction Lagrangian 

C mi = -i9 ■ = -ga a^-d^n. (19) 

This leads to the so-called 7r — a\ mixing, which is usually eliminated by a redefinition of 
the a\ field, ai^new = ai^oid + gco/m^ d^n , followed by a wave- function renormalization 
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Figure 1: The pirn vertex function of the gauged linear a model after symmetry breaking. 



of the pion field, 7f new = Tf \d/y/Z, with Z = m 2 ai /m 2 p . As a consequence, the physical 
pion mass is increased by y/Z, 

= ^-m^o , (20) 
m p 



whereas the pion decay constant is reduced by l/y/Z, 

U = ^<r Q . (21) 

m ai 

However, we will choose a different method to treat the it — a± mixing, namely summing 
the self-energy diagrams generated by the it — a\ vertex to all orders, which gives the same 
results at tree level, but is more convenient if loop corrections are taken into account. In 
the next section this will be discussed in more detail. 

The properties of the p and the a\ are dominated by the decay modes p — > tttt and 
ai — > irp, respectively. The corresponding vertex functions within the present model are 
shown in Figs. [I] and |2[ If the i\ — a\ mixing is not eliminated, there are four contributions 
to the pirn vertex. The corresponding vertex factors read: 

{rfiy abc = -ge abc (2k» + qn, (22) 
(r^ c) )l = $ i ^ c (2F + ^), (23) 

$%)L=^ ( 24 ) 



The coupling constant g can be determined from Eqs. QTS[) and ( |2"T| ) and the empirical 
values of f n , m p and m ai . If we take e.g. m p = 770 MeV, m ai ~ \f2m p = 1090 MeV , 



/■ k + q,v,c 




(a) (b) 

Figure 2: The aiirp vertex function of the gauged linear a model after symmetry breaking. 
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and f n = 93 MeV, we find g = 5.9 . This is more or less the value which is needed 
to obtain a good description of the p-meson width and the electromagnetic pion form 
factor in models without ai-degrees of freedom, where the pirn vertex function consists of 
diagram (a) only (see e.g. Ref. |§). In the present model, the second contribution, (b+c), 
has the same structure as (a) and leads to a reduction of the pnir coupling constant by a 
factor 1 — g 2 a\jrp? ax = rp? p jrr? ai = 1/Z. However, this is exactly cancelled by the factors 



Z assigned to each external pion line if one calculates the decay p — > nn. Thus, the 
contributions (a) to (c) would give a quite good description of the decay p — > 7nr with 
the value g = 5.9 determined above. Unfortunately there is still contribution (d), which 
generates a strong g 2 -dependence of the vertex function. This does not only result in a 
too small width for p — > -kit, but also leads to a wrong shape of the pion electromagnetic 
form factor or the ir — n scattering phase shifts in the J = 1,1 = 1 channel as a function 
of q 2 , which can perfectly be reproduced without momentum dependence of the vertex. 

A similar effect exists for the decay a,\ — > irp. Due to ir — a x mixing there are two 
contributions to the a^np vertex function shown in Fig. |2|: 

KUI = ^V*^. (25) 
(r£X c = -i^- s abc (Pr + fcV + - (k 2 + 2 k ■ q) sT) ■ (26) 

m a 1 



Again, the term arising from the non-abelian terms in the field strength tensors Y, 



Cyyy = -|tr^[^,n (27) 

causes an additional momentum dependence of the vertex. As a consequence, the ampli- 
tude for ai — >• pir as a function of the invariant mass of the at has a zero at q 2 = m 2 ai +m 2 . 
This lies well below the r mass, i.e. within the range where the a\ can experimentally be 
observed under clean conditions. The experimental data, however, show no minimum in 



the Tip invariant mass spectrum |21, |2 



2.2 Restriction to global chiral symmetry 

The problems discussed in the end of the previous section force us to modify the La- 
grangian, Eq. ([15]). To that end we relax the constraints imposed by the requirement 
that the Lagrangian (except for the mass term C mo ) should be invariant under a local 
SU{2)l x SU(2)fi transformation and only insist on a global chiral symmetry. In fact, 
this has already been done in Ref. |l9| , but in a less radical way. 

Two consequences of the gauge principle are incorporated in the Lagrangian (|T5j): 
First, the same coupling constant g appears in different interaction terms, and second, 
some chirally symmetric interaction terms are not allowed. For example, the most general 
interaction of two $ fields with two Y fields which is chirally symmetric and has dimension 
< 4 looks as follows: 

C^yy = ~ y Y^ ■ Y"$ + ^ $ • $ tr Y^ . (28) 

The gauge principle would force us to choose hi = g 2 and h 2 = 0, whereas in Ref. |19[, hi 
and /12 were chosen as independent parameters hi = g 2 (1 + 2c) and h 2 = g 2 (b — c) (with 
b and c as defined in Ref. [0). 
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As we have shown in the previous section, the problems with the phenomenology 
of p and ai mesons arise mainly from the vector meson self interaction. However, if 
we are interested only in global chiral symmetry, we can adjust the coupling constants 
independently: 



Cyyy + Cyyyy = ~Y tT d * Y » [ y "> Y "] + X tr Y ^ [ y "> Yl 



(29) 



We are allowed to choose the coupling constants g' and g" smaller than g or even to put 
them to zero. In fact, we perform our calculations with arbitrary values of these coupling 
constants, and find that the best description of the axial- vector correlator can be achieved 
with g' = 0, while our results are insensitive to g" . For simplicity, we will put both, g' 
and g", to zero throughout this paper. 

In addition to the g" term, one could think of other chirally symmetric YYYY- 
interaction terms of dimension 4. However, they are not important for our purposes, 
either. Higher- dimension terms, which are frequently used in the literature [|TI| [T?|, [19| 



are much more complicated and result in momentum dependent vertices. As we will see, 
the experimental data on p and ai mesons can very well be explained without such terms. 
This indicates that in nature the momentum dependence of the vertices is small. 

Besides the interaction terms, we also modify the kinetic term for the vector mesons 
contained in Eq. ([13|) . With the original kinetic term, the longitudinal part of the vector 
meson propagator does not go to zero for k 2 — > oo, causing serious problems if one tries 
to calculate loop diagrams with vector meson lines. To avoid these problems, we take the 
Stiickelberg Lagrangian 



Adnr = -\ tr [d,Y v - d v Y^Y v - - | tr {d^f . (30) 

From this one retrieves the kinetic term of the original Lagrangian (O) by taking the 
limit £ — > . The p meson propagator now has the form 

k^k" _ nv k^k v 
TnP 1 Trip" 

1 k z — m z p k — rrip/E, 

This form has the advantage that for £ > the propagator behaves as 1/k 2 for k 2 — > oo, 
which is not true for £ = 0. On the other hand, we have introduced unphysical scalar 
particles with mass m p j \/f . This can be interpreted as some form of regularization of the 
contributions from the longitudinal polarization states of the vector mesons. Therefore 
the parameter £ must be chosen small enough, such that all thresholds involving these 
unphysical particles lie above the energy range we are interested inF]. 
Combining everything, we obtain our full Lagrangian: 



1 u 2 h h 

- ■ <9 M $ - — $ ■ $ - ig r M $ • 6^$ - — r M $ • + ^ 

A 2 



($ . $)2 + ca _ 1 tr ( d Yv - d v Y„)( d»Y v - d v Y») - | tr {dpY^ f 
4 o 4 

2 

+ ^tr^. (32) 



^^In a true gauge theory the final results are independent of £, because in this case the contributions 
of the unphysical scalar particles are cancelled by the ghost contributions. 
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7T 7T 7T 0-1 TT 

(a) k = + •UZQQ^ 

(b) k, // vAfififlfiflfijb ' = ^QQQQQQQflv + vfi&L»- ^QflQ> 

TT ai 

(c) A: -> A« 

Figure 3: Propagators with 7r— ai mixing: (a) the pion propagator iG 7r (k), (b) the ai propagator 
iGai(k), and (c) the mixed propagator iG% ai (k). 



Again, as in the previous section, we must shift the a field, and Eqs. fll6|) and fll7|) 
remain valid. But instead of Eq. flUf ) we find the following expressions for the meson 
masses: 

2 2 , \2 2 2 2 , o\ 2 2 2 2,7,2 2 2 , tr. , i, \ 2 

m n0 = /i + A a Q , m CT = /i + 3A a , m p = m + h 2 a , m ai = m Q + [hi + h 2 ) a . 

(33) 

This has the interesting consequence that with a special choice of parameters (mo = 0, 
h 2 = m 2 p ja^) it is now possible to generate the vector meson masses from chiral symmetry 
breaking alone. In this case the masses of both p and a\ meson would drop in a nuclear 
medium where chiral symmetry is (partially) restored ("Brown-Rho scaling" [[|). 

The 7i — ai mixing Lagrangian C nai is still given by Eq. flTD|), but because of the 



Stiickelberg (£-) term in the Lagrangian it cannot be eliminated by a simple field redefini- 
tion. Instead, self-energy contributions from C nai must be iterated to all orders as shown 
in Fig. 0(a). Then the pion propagator can be written as 

G -W= ... 1 ?a 'Eict' < 34 > 



]A - m 2 4- g2fT o fc2 ^ k 2 -m 
ft "SrO + Cfe 2 -™^ »=1 



TT / 



The first pole at m„-i = m n corresponds to the physical pion, whereas the second one at 
m n2 corresponds to the unphysical longitudinally polarized at particle. In principle m^i 
and m 7T 2 are functions of m n o > but it is more convenient to quote the relations for these 
masses as functions of the physical pion mass: 

2 2 2 m lS m l, -3 2 °l - W) 2 m >a! ~ /nr\ 

77^ = 77^, m l2 = — -= — ^ , m w0 = — — . (35) 



The Ziri are given by 



^(m^-^) " m 2 ax -iml 



Z n l = —, \ 2~T ' ^tt2 = 1 — ^ttI • (36) 



Om\ 2 - m 



7T 1 / 



Here Z^i corresponds to the Z factor introduced in the previous section. In the limit 
^ — > 0, /ii — ^ (? 2 and h 2 — > we retrieve the relations for m n , m n0 and Z in the gauged 
linear sigma model. 

The a± and the mixed (tt — a±) propagator can be obtained from the pion propagator 
as shown in Figs. ||](b) and (c). For the a\ propagator we find 

GZ{k) = ^ + WW T^V > (37) 

1 k z — mi, k — mi: 

1 i — \ 
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with 



2 2 

ml , — mi 



7 - 'ZJll n0 7 - 7 Cv 

A ai 1 — 7 — h~l 2 2~T ' a i 2 ~~ 2 a i 1 ' V 



^2/2 2 \ ' 'CU z 

^m 2 1 (m 2 w2 -m 2 nl ) m ; 



and for the mixed propagator 
with 



2 7 

A; 2 — m 2 
i=i 7rl 



-^Trai 1 — — 77 a 2~T ' Z nai 2 — — Z^ai 1 • (40) 



m 



Since we have switched off the 3-vector interaction, for the pnir vertex functions only 
the diagrams (a) to (c) of Fig. [I] are left. If the external pions are on-shell, i.e. k 2 = 
(k + q) 2 = m\ , the total vertex function reads, 

( r ^)l = ~{g- m2 ^L ) e abc (2 P + <f) . (41) 
This can be used to define an effective pnn coupling constant 



e« _MLU _ K~eO«-^) 

9pmr — [9 o C 2~) 2 C 2A2 2~~2 T " ^ J 



2.3 Electroweak interactions 

Until now we have completely neglected the electroweak interaction. However, the cleanest 
experimental information about p and a\ mesons can be obtained from electromagnetic 
reactions, e + e~ — > p, and from the weak decay of the r lepton, t u T p or t —>■ u T a\. 
Another reason why we want to include the electromagnetic interactions is of course our 
interest in dilepton production in heavy ion collisions. 

At the quark level there is no ambiguity how the electromagnetic field and the W 
bosons must be introduced: The ordinary derivative of the quark fields is replaced by a 
covariant derivative. If we consider only u and d quarks (ip = (u, d)), this looks as follows 



(see e.g. [gl): 



n / fa ■ a r3 ^ e cos #c 1 - 7s n t 2 \ \ 

m = - "4* 2 - — r 1 " 2 + ^ 2 J + • • -r ■ (43) 

Here we have omitted the isoscalar part of the electromagnetic interaction and the Z 
boson coupling, since they are irrelevant for our purposes. 

This can be generalized to the (a, 7?) meson sector. The isospin operator r/2 cor- 
responds to the T operator, whereas 75T/2 can be identified with T 5 . Therefore all 
derivatives <9 M $ in the Lagrangian fl32|) have to be replaced by 

ie cos Qr 



Dfi = (d, - ieA„T 3 - {W lf jf + W 2 »T 2 L ))$ 

V Sill uw ' 

9 v (T ) + e ( ~ ° ) + eC ° s6c ( ^ ' 1 

<9 M 7?y I X x 7?/ 2 sin 9 W \W, A x 7? + a 



(44) 
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In the last line we have used the abbreviations 




4, = | | , = | W 2fl | . (45) 

To find the operators in the vector-meson space corresponding to the T and T 5 operators 
in the (a, 7r)-meson space, we can expand Eq. flltf) for the case d^U = to linear order 
in a and [3. Then it becomes clear that the derivatives of the vector meson fields, d^Y v , 
must be changed into 

DJT V = djf v - te [T 3 , K] - WCOs 6c (W ltM [T\ L , Y v ] + W 2fl [T 2 L , Y v \) 

Sill Oy/ 

=:D hl p u -f + D^a lu -f\ (46) 

with 

j— . ^ „ ^ t* ^ ecos^c - _ v 

= dpPv -\- e A^ x p v -\- , x (p„ - ai „) , 

z sin (7^/ 

— * 6 COS 

D^a lu = d^a X v + e A^ x a lu + . W M x (a lu - p v ) . (47) 

l sin t7]4/ 

One should mention that the replacements described above must be made before the 
redefinition of the a field, Eq. (fL6D - When we now perform the shift of the a field, W — tt 
and W — a± couplings emerge. As a consequence there are two contributions to the pion 
decay, a direct one and a contribution with an intermediate longitudinal a±: 

f , = JZlU - -r&j-r) = ;°'<-^-^» . ,48) 



Obviously we did not assume anything about vector meson dominance when we con- 
structed the coupling of the A and W fields to the mesons, and in principle the vector 
mesons were not treated differently from the scalar mesons. Hence, the A and W fields 
simply couple to the Noether currents. In addition there are seagull terms involving two 
A or W fields. At this stage there is no direct 7 — p or W — p coupling in our model. 



However, as we will see in Sect. fO, the 7 — p and W — p coupling appears naturally at 



the one-loop level, where we will need a counterterm of the form 

C, Y + C WY = -{ tr (d,Y v - d v Y,) (eA u T 3 + (W^T? + W»T 2 L )) . (49) 

4 V Sin U\y / 

The 7 — p part of this term has exactly the form of the 7 — p coupling of Kroll, Lee and 



Zumino |25] 



£yp = - Y (<9 ^ ~ d "P*^ AV ~ dVA ") > ( 5 °) 



the only difference being that in Ref. [25 the coupling constant / was fixed to 1/g, whereas 
it is a free parameter in our model. 
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3 Meson propagators in one-loop approximation 



Our aim is the chirally symmetric description of the vector isovector and axial-vector 
isovector correlation functions. This implies that we need a chirally symmetric description 
of p and d\ mesons. By construction, our Lagrangian is chirally symmetric, but most 
approximations beyond tree-level destroy the symmetry. For example, in the linear a 
model without vector mesons the mean field (Hartree-Fock) approximation violates the 
Goldstone theorem, which can be recovered only if one includes also the RPA corrections 
p6| . Then, of course, the self-consistency is lost, i.e. the final RPA mesons are different 
from the mean-field mesons propagating in the loops. In fact, a more general investigation 
of this problem shows that there seems to be no self-consistent approximation scheme 



which respects the symmetry at the correlator level ||27|| . 

A simple symmetry conserving approximation scheme is given by a loop expansion, 
which formally corresponds to an expansion in powers of ft. Here we will calculate the 
meson self energies to one-loop order. This approximation has the following properties: 

1. The pion remains massless in the chiral limit (Goldstone theorem). 

2. All differences of chiral partners (it — a and p — a\) disappear when the chiral 
symmetry is restored (<t = 0). 

3. The dominant meson decay modes, such as p — > tttt or a\ — > np, are contained in 
the scheme. Of course, due to the lack of self-consistency the final decay of the p 
meson in a% — > Tip — > tttttt is not included. 

Of course, in a strongly interacting theory the loop expansion does not converge. 
Therefore the model is not defined by its Lagrangian alone, but only by the Lagrangian to- 
gether with the approximation scheme used. Hence the propagators and vertices obtained 
from the Lagrangian must not be interpreted as elementary, but rather as parametriza- 
tions of the full propagators and vertices, containing all non-perturbative effects which 
are not explicitly taken into account in the approximation scheme. 



3.1 Subtraction of divergences 

If one calculates diagrams with loops, one always finds divergent expressions which have 
to be regularized in some way. The above mentioned fact, that the vertices and prop- 
agators obtained from our Lagrangian are not elementary, justifies the introduction of 
form factors, which are frequently used in the literature in order to cut off high momenta 
0, |3C§, rendering the loop integrals finite. These form factors are independent pa- 
rameters and cannot be computed from the Lagrangian. Unfortunately we cannot simply 
attach form factors to the vertices obtained from our Lagrangian, since this would destroy 
the symmetry properties of our model. 

Other regularization prescriptions, which can be found in the literature, are the Pauli- 
Villars scheme with finite regulator mass [|5] or subtracted dispersion relations. Using these 
prescriptions one can sometimes avoid the problems arising from form factors, although 
they give very similar results. As an example consider the self-energy Ti p _> nn , describing 
the decay of the p meson into two pions, with a three dimensional form factor v(k) 
attached to each pirn vertex, which cuts off pion momenta k (measured in the p-meson 
rest frame) at some value |/c| max . For small p-meson energies, q 2 < 4(m^ + fc^), the 
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imaginary part of Sp-^ is not modified by the form factor. For the real part the situation 
is more complicated. Since it is related to the imaginary part via a dispersion integral, it 
receives already for low q 2 contributions from the imaginary part at high q 2 . Therefore it 
depends strongly on |fc| max - However, using subtracted dispersion relations, one can easily 
convince oneself that for low energies the influence of the form factor on the real part can 
approximately be absorbed in the subtraction constants, while the remaining dispersion 
integral is almost independent of the form factor. 

In this sense, subtracted dispersion relations can be interpreted as an approximate 
way of attaching form factors to the vertices. As the subtraction constants depend on the 
high-energy behavior of the form factors, they cannot be calculated from the Lagrangian, 
but they are free parameters. However, there are some contraints which follow from chiral 
symmetry. In order to preserve the symmetry, we must make equal subtractions from the 
self energies of chiral partners. Obviously this is mathematically ill-defined, since the 
subtractions are infinite. Therefore we need a more systematic method to ensure the 
consistency of the subtractions. 

The subtractions can be regarded as (infinite) counterterms added to the original La- 
grangian. Thus the subtractions are consistent with chiral symmetry if the counterterms 



are chirally invariant. The values of the counterterms are determined as follows [p8 [. First 
we calculate the self-energy diagrams with dimensional regularization, i.e. we calculate 
the loop integrals in d space-time dimensions, 



d A k fd d k /A 2 v, 

(51) 



- 



2 



(2tt) 4 J (2ir) d \AttJ 

If we expand the result in the vicinity of d = 4, the divergences appear as poles oc 1/(4— d). 
The residues of these poles are mathematically well-defined, and the self energy can be 
written in the form 

e = s+e ~(iVD- < 52 > 

In this way we have split the self energy into a finite part E, and an infinite part. (7 = 
0.577. . . = Euler's constant is finite, but it is convenient to keep the infinite 1/(4 — d) 
term and the 7/2 term together.) The counterterms are also split into finite and infinite 
parts, e.g. 5fi 2 = 5fi 2 + 5^(1/(4 — d) — 7/2). The infinite parts of the counterterms must 
be adjusted such that E ro vanishes, while the finite parts are treated as free parameters. 
The crucial point is that each counterterm may be adjusted only once and not separately 
for different self energies. 

The factor (A 2 /(4vr))( 4 - d )/ 2 in Eq. (0) ensures that the dimension of the result does 
not depend on d. The finite parts of the integrals, of course, depend on the choice of A. 
However, the physical results are independent of A, if the finite parts of the counterterms 
are appropriately readjusted after a change of A. For example, if A is changed to A', 5fi 2 
must be changed to 5jj? + dji 2 ^ ln(A'/A). 



3.2 The tadpole correction 

At tree level, no a line can disappear into the vacuum, because all terms linear in a 
are eliminated by the shift of the a field, Eq. (^6|), as a consequence of Eq. (|i~7|), which 
can be illustrated diagrammatically as shown in Fig. ^(a). This is no longer true at the 
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Figure 4: (a) Diagrammatic illustration of Eq. (]i~7|), which shows that there are no tadpole 
graphs at tree level, (b) Tadpole graphs T at one-loop order. The meaning of thick and thin 
lines is the same as in Fig. |3[ 



one-loop level, where one has the additional tadpole graphs shown in Fig. |4](b). The last 
two contributions are generated by the chirally symmetric counterterms 

SCp + 5C X2 = — ^$ ■ $ - — ($ ■ $) 2 . (53) 

Calculating the diagrams displayed in Fig. ^(b) and splitting the result into finite and 
infinite parts, T = T + Too(l/(4 — d) — 7/2), we find for the infinite part 

^ = 3g, (- (fti + 2WK _ ffV + (3 , 8i + 6ft2)m? + 2A2(i2 

+ 3^ + g y + ft j _ 2 g * {hl + h2 ) + + + + iX , 

-iao5fil -ia^5X 2 00 . (54) 

From this one can immediately see that T can be made finite (i.e. = 0) by a proper 
choice of Sf/^ and SX 2 ^ . Furthermore, this choice does not depend on c or a , i.e. on the 
explicit or spontaneous symmetry breaking. In the following sections we will see that this 
is a general property of all counterterms, as shown by Lee [25] for the linear a model 



without vector mesons and by Ram Mohan |28| for the linear a model with nucleons at 



finite temperature and density. In this context we should remark that, in the limit £—►(), 
it is impossible to find suitable counterterms independent of a . For example, the third 
contribution in Fig. |](b) can be written as 

-3k ao g^Jj^GFik) = 3h 2 a f±L + _L_) 

= 3 X>2 J^y ( (jfc2 _ m 2 )n + (efc 2 _ m 2 )n ) • ( 55 ) 

For £ > the n = 1 and n = 2 terms are quadratically and logarithmically divergent, 
respectively, while all terms with n > 2 are finite. This is the reason why two counterterms 
independent of <r are sufficient to cancel the divergences. In contrast, for £ = 0, all 
terms are quartically divergent, i.e. one needs either an infinite number of counterterms 
independent of a , or counterterms which depend on <j Q . 
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Figure 5: One-loop contributions to the p-meson self energy, — iS p . 



3.3 The p-meson self energy 

In Fig. [5] all one-loop diagrams for the p-meson self energy are displayed. The last three 
diagrams correspond to the contributions of the three counterterms which are needed to 
cancel the divergences: 

6C m 2 + 5C Zy + 5C h2 = — ± tr + -f tr (d M Y„ - d u Y,) ( W - 8TY») 
u 4 o 



Sh 



+ — ifc.fctrY^y (56) 



4 



Again, the self energy is split into finite and infinite parts, = E p + E^(l/ (4 — d) — 
7/2), and the infinite parts of the counterterms are chosen such that E^ vanishes^ . 

Since in our model the p meson is not a gauge boson, the p-meson self energy is not 
transverse. Instead, it can be decomposed into transverse and longitudinal components, 

= K(i 2 ) <cf - 9" u ) + K^^-f- (58) 

As we will see, the longitudinal part does not contribute to the vector current correlator 
as a consequence of current conservation. 



3.4 The pion and 01-meson self energies 

Because of the 7r — a\ mixing, the situation for the pion is a bit more complicated than 
for the p meson. It is necessary to distinguish between the one-particle irreducible (1PI) 
self energy E^ , which cannot be cut into two parts by cutting any single line, and the 
proper self energy , which cannot be cut into two parts by cutting a single pion line, 
but which can possibly be cut into two parts by cutting a single ai-meson line (see e.g. 
the last diagram in Fig. 0(a)). 

For the a\ meson the situation is similar. Here we have to distinguish between the 1PI 
self energy E aifll , which cannot be cut into two parts by cutting any single line, and the 

2 For non-vanishing vector-meson self interaction, i.e. g' 7^ or g" 7^ 0, there are more diagrams than 
shown in Fig. and one additional counterterm is needed: 

<5£ e = -|tr(5^) 2 . (57) 
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Figure 6: One-loop contributions to the 1PI pion self energy, — iY> T , 



proper self energy £ ai , which cannot be cut into two parts by cutting a single ai-meson 
line, but which can possibly be cut into two parts by cutting a single pion line. 

In Fig. |B| we display the one-loop contributions to the 1PI pion self energy H n7T . The 
last three diagrams are generated by the counterterms given in Eq. (|53|) and the additional 
wave-function renormalization counterterm, 



2 M 



One can explicitly show that this self energy fulfills the Goldstone theorem, i.e. 

S^(0)=0 for c = 0. 



(59) 



(60) 



In particular, it is interesting to see that in the chiral limit the contributions from the S/j, 2 
and SX 2 counterterms (14 th and 16 th diagram) exactly cancel the counterterm contribu- 
tions contained in the tadpole correction (13 th diagram). 

The one-loop diagrams for the 1PI ax-meson self energy are displayed in Fig. 

[7|. The diagrams 9, 10 and 12 are generated by the counterterms given in Eq. (|5"rf ). In 
addition we need one more counterterm, namely 



5h 
2 M 



(61) 



which gives rise to the 11 th diagram. In analogy to the p-meson self energy (cf. Eq. (jj8|)), 
the ai-meson self energy can be decomposed into transverse and longitudinal components, 
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Figure 7: One-loop contributions to the 1PI ai-meson self energy, — zS aiai . 
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Figure 8: One-loop diagrams for the 1PI mixed ir — a\ self energy, — iS^ ai . 



E* and E^ . The transverse part is not modified by it — a± mixing, and therefore the 
proper transverse ai-meson self energy is equal to the 1PI transverse self energy, 

KM 2 ) = KaM 2 ) • (62) 

To obtain the proper pion self energy, we have to consider also the 1PI mixed self 
energy E£ a , which is a correction to the it — ai vertex given by the Lagrangian (p!9[). The 
corresponding one-loop diagrams are shown in Fig. [8[ The last diagram comes from the 
counterterm 

5C g = -i5g ■ . (63) 
With this mixed self energy, the corrected 7r — a\ mixing vertex is given by 

(IW?))^ = -5 ab (g* q» + i^ ai (q)) =■ -<f<U F„ ai {q 2 ) , (64) 

where q is the incoming momentum of the pion. The proper pion self energy contains the 
tree-level self energy shown in Fig. |3](a) and the loop corrections, i.e. 

S.(g 2 ) = E^(g 2 ) + - q2 [ F : aiiq y\ 9 , . (65) 



m 2 

ai 



The corrected pion propagator is given by G4 = l/(<? 2 — m 1v ~ ^n)- ft * s obvious that 
T, n fulfills the Goldstone theorem if E T7r does, because the contributions due to n — a\ 
mixing are proportional to q 2 . 

3.5 The (T-meson self energy 

The self energy of the ex- meson, E CT , is given by the diagrams shown in Fig. ||]. To make 
it finite, the counterterms Sfi 2 , SX 2 and 5Z$, introduced already in the previous sections, 
are sufficient. 

4 Vector and axial-vector correlation functions 

Let jn and be the hadronic vector isovector and axial-vector isovector currents, re- 
spectively. In terms of quark fields ip = (u,d) these currents read = and 
Ju = V'7m75^V'/2- The current-current correlation functions are defined by 

K\M = * Jd 4 xe^(0\T(f a (x)f b (0))\0) , (66) 
n^(?) = i |rf 4 xe^(0|T(j^(x)jr(0))|0) . (67) 
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Figure 9: One-loop contributions to the u-meson self energy, — iS c 



Because of isospin symmetry we have ^y/ Aab 



Tly, A 6 a b . The correlation functions can 



be split into longitudinal and transverse components as described for the p- and ai-meson 
self energies: 



(6c 



Since the vector current is conserved, its longitudinal correlator vanishes, i.e. H v = . As 
usual, inserting a complete set of energy and momentum eigenstates into Eqs. (|66"D and 
(|57|), one can derive spectral representations for the correlation functions. Following the 
definitions of Ref. p2| |, the corresponding spectral functions v i, a± and are given by 



v±(q 2 



An 
q 2 



ImlT^g 2 ), 0l ( 



Air Att 

— Imn* A (g 2 ), a (q 2 ) = —lmU l A (q 2 ) . (69) 



To order e 2 , the current-current correlation functions are directly related to the photon 
and VT-boson self energies: 



yfiu 



y[iv 



—e 2 U^ + isoscalar contributions 



e cos 9 c 
2 sin 9w 



;nr + iFf ) 



Experimentally, Ily can be measured in the reaction e + + e — > hadrons. Up to q 1 
both ily and 11^ can be measured in the reaction r — ► z/ T + hadrons. 



(70) 
(71) 

9 



4.1 Photon self energy and vector correlator 

In this section we will calculate the spectral function V\ of the vector correlator, which is 
proportional to the imaginary part of the (isovector) photon self energy £ 7 . Similar to 
the pion and ai-meson self energies, the proper photon self energy £ 7 has a one-particle 
irreducible (1PI) contribution, £ 77 , and contributions which can be cut into two parts by 
cutting a single line, which in this case must be a p-meson line. 



In Fig. [10| we display the one-loop diagrams for the 1PI photon self energy, S^ 7 . Here 
we do not introduce counterterms to cancel the divergences, because we are only interested 
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Figure 10: One-loop diagrams for the 1PI photon self energy, —iT,^ u 



77- 



in the imaginary part, which is finite. The last three diagrams are purely real and do not 
contribute to the spectral function v\, but taking them into account one can explicitly 
show that within the dimensional regularization scheme the self energy is gauge invariant, 
i.e. 

q^ ( (q) = and E£(0) = 0. (72) 

For the remaining self energy contributions we must compute the 7 — p transition 
vertex. The corresponding mixed self energy, E^(g), is given by the diagrams shown in 



Fig. |TT| Again, the conditions for gauge invariance, 

q^(q)=0 and E£(0)=0, ( 73 ) 

are fulfilled. This is the reason why only one counterterm is needed to cancel the remaining 
divergence, namely a counterterm of the same form as the direct gauge-boson — vector- 
meson coupling, Eq. ([49|), 

S£f = - S 4 tr {d,Y v - d v Y,) d» (eA u T 3 + (W?T{ + W»T 2 L )) . (74) 

4 V sin t>w / 

In terms of the mixed self energy, the one-loop corrected 7 — p vertex function can be 
written as 

(r ip (q))7 = iS 3 a(ef(q"q v - - =: iefoCW - 9V) F lp (q 2 ) . (75) 

Here a is the isospin index of the p meson. 
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Figure 11: One-loop diagrams for the mixed 7 — p self energy, — iT, 
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Figure 12: One-loop diagrams for the 1PI axial-vector VK-boson self energy, — (seagull 
and tadpole graphs omitted). 



Now we can write down the total photon self energy. Because of Eqs. ([72]) and (|73|) 
its longitudinal component vanishes. The transverse component reads 



<76) 



Finally we obtain the vector spectral function V\(q 2 ) by taking the imaginary part of S* 
see Eqs. © and ©. 

4.2 Transverse part of the axial- vector correlator 

Similar to the vector correlator, the axial-vector correlator has one-particle irreducible 
contributions, and contributions which can be cut into two parts by cutting a single 
line, in this case either an ai-meson or a pion line. Of course, the single-pion states can 
contribute only to the longitudinal component of the correlator, which will be postponed 
to the next section. 

Due to the parity violation of weak interactions, the W-boson self energy has vector 
and axial-vector contributions, Ey^ = Sj^ + . In this section we are looking only at 
the axial- vector contribution, S^. Its one-particle irreducible (1PI) part, X^Vjy, is given 
by the diagrams shown in Fig. |T2| and some real seagull and tadpole contributions, which 
we have omitted since we are only interested in the imaginary part. As in the previous 
section, the imaginary part is finite. Since the axial current is not conserved, * s no ^ 

transverse and can be split into transverse and longitudinal components, and ^ww- 

In Fig. 13 the mixed W — a\ self energy, , is shown. The missing current conser- 
vation does not only lead to a non- vanishing longitudinal component, T, l Wa , it also avoids 
the cancellation of divergences due to the condition S|^ a (0) = 0, which would be valid 
if the current was conserved. Hence, the single 6f counterterm given in Eq. ([73]) is not 
sufficient to render S^ ai finite. However, the minimal substitution, Eq. (fHj), must also 
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Figure 13: One-loop diagrams for the mixed W — oi self energy, — i£^ al . 
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be applied to the Sg counterterm introduced in Eq. (pB|). This finally leads to the last 
diagram in Fig. 13, and indeed, with this counterterm the remaining divergence in S^ ai 



is cancelled. The total 1-loop corrected W — a\ vertex can now be written as 

ecosfl c ((q»q v \ t 2 ^ , 2 \ 

Here a is the index of the W field (1 or 2) and b is the isospin index of the a\ meson. 
The spectral function of the transverse axial correlator, ai(q 2 ), is now obtained from 

T At,^_ T At r „2x | ( ecos6 c \ 2 (^ Ql (g 2 )) 2 / ? M 

V2smf%/ g - - 2^ 

together with Eqs. and ([7T|). 

4.3 Longitudinal part of the axial- vector correlator 

For the longitudinal part of the axial- vector correlator, we must also consider contributions 
with single-pion intermediate states. To that end we need the W — it transition vertex, 
which again consists of a one-particle irreducible (1PI) contribution and a contribution 
which can be cut into two parts by cutting a single ai-meson line. The 1PI diagrams for 



the mixed W — ir self energy, are displayed in Fig. [14]. The last diagram is generated 
by the counterterm one obtains by performing the minimal substitution given by Eq. (13) 
in the wave-function renormalization counterterm, Eq. floTf). The 1PI one-loop corrected 
W — 7T vertex can now be written as 



= + = : Sf SM2) • (79) 

with q being the incoming momentum of the W boson. Then the full W — ir vertex reads 
FwM)* - -J—eZ* 5 ab[F W7r (q ) - —————— 



2sm9 w V rn^-iq 2 -T} aiai {q 2 ). 

q^ ab F^ n (q 2 ) . (80) 



ecos0 c r tot/2N 



2 sin 9\y 



21 



The longitudinal iy-boson self energy, which determines the spectral function ao{q 2 ) ac- 
cording to Eqs. (]69|) and (]71~D , is given by 



T Al (n 2,_yAl ,»s | f eCOsOc Yf {F l Wai {q 2 )) 2 , ^WW)) 

^C? j - W? J + Usin^J Kml - - Si iai (g 2 ) + q* - m 2 - 



*1) 



The full W — tt vertex can also be used to calculate the one-loop corrected pion decay 
constant, 



= FfcirnP*) y/ Z?\ (82) 

with mi and Z-tK defined as 



nrj - =ml + E^m^ 2 ), (83) 



,(1)2 



5 Results 



5.1 Determination of model parameters 

Our model has two groups of parameters, the bare parameters of the Lagrangian, 

/i 2 , A 2 , c, f, mo, fci, /i 2 , /, 

and the finite parts of the counterterms, 



8Z§ , 5/j 2 , 5X 2 , bZ Y , 5ml, Sg, Shi , 5/ • 

However, not all of these parameters are independent of each other. Our final results 
do not separately depend on / and 5f, but only on the combination (/ + 5f), and the 
counterterms 5ml and 5h 2 always appear in the combination (5ml + 5h 2 al). Thus we 
have altogether 16 independent parameters. 

In practice it is more convenient to employ Eqs. (|TT|), (|33"D, (|35|), (f42|), and ( f4~8"|) and to 
express the bare parameters p 2 , A 2 , c, g, /ii , and h 2 through the more physical quantities 

m n , m CT , m p , m ttl , , and ^ . 

For the parameter fixing one should keep in mind that our approximation is not self 
consistent. In particular, unstable particles calculated in one-loop approximation decay 
into tree-level particles. Hence one should not only fit the one-loop properties to the 
empirical values, but also the tree-level masses, in order to get the proper phenomenology. 
For instance, the decay a\ — > txtxtx is modeled by the decay modes a x — > Tip (dominant) 
and a± — > na (small), with stable p and o mesons. As we shall see, for the Tip decay 
mode, this works rather well for invariant masses above m„ + m p , if both, the pion and 
the p meson, at tree level have already their physical masses, i.e. m n = = 140 MeV 
and m p = m p . As a definition for the p-meson mass we choose the position where 
the p-meson spectral function has its maximum, m^ = 767 MeV. 
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Besides the pion mass we also fit the pion decay constant to its empirical value, 
f w = 93 MeV. Again we demand that the tree-level value is equal to the corrected one, 
f n = fjp. Among others, this will be important for later applications of the model at 
non-vanishing temperatures, where the effects will be dominated by thermally excited 
tree-level pions. 

Similarly, we demand that m ai and m a at tree-level coincide with the corrected a\- 
and cr-meson masses, i.e. m^} = m ai and rria = rn a . The corrected masses are again 
defined as the positions where the spectral functions have their maxima. As we will discuss 
below, one should not blindly fit the ai-meson mass to any value given in the literature. 
Instead, we will use the ai-meson mass to fit the peak position of the axial vector spectral 
function ai(q 2 ). For the cr-meson, which is not a well established resonance, we will take 
two different values, namely 600 MeV and 800 MeV. 

The width of the p meson depends on two parameters, and 5Zy ■ The latter has 
also influence on the a\ meson width. Hence, both parameters can be fixed by fitting the 
widths of the vector and axial- vector spectral functions, vi(q 2 ) and a\(q 2 ). 

The width of the vector spectral function v\(q 2 ) being fixed, its height can be adjusted 
with the parameter (f + Sf). The height of the axial vector spectral function ai(q 2 ) is then 
adjusted by changing 5g, which, in practice, is quite complicated, since this parameter 
changes also many other quantities because of tt — a\ mixing. 

The parameter £ acts as a cutoff for the momenta of longitudinally polarized vector 
mesons. If we set arbitrarily £ = 0.1, the masses of the longitudinally polarized vector 
mesons are > 2 GeV, which lies still in the range of hadronic scales, but is high enough 
to avoid unphysical thresholds in the energy range of interest. Since the masses of the 
longitudinally polarized vector mesons behave as 1/ we can obtain reasonable fits also 
for smaller values of £, e.g. for £ = 0.025. 

The remaining undetermined parameter is mo, the vector- meson mass in the restored 
phase. This value is partially constrained by the upper limit of the a\ — > no branching 
ratio. The relation between m and the a\no vertex function T aina is the following. T ai7T(7 
has two contributions, corresponding to the diagrams (a) and (b) in Fig. |2|, with the 
outgoing p-meson lines replaced by cr-meson lines. The sum of both terms yields 

(rK>): t = -^(^^2F + ,-). (85) 

For transversely polarized a\ mesons, the c/ 1 term does not contribute. Thus, for £m 2 <C 
m 2 av the partial decay width a± — > no is proportional to m\ . However, even for m = the 
axial- vector spectral function a\{q 2 ) has small 7rcr contributions from the direct W — n — o 
coupling (first four diagrams in Fig. |l^). 

It turns out that we can achieve reasonable fits for mo < 400 MeV. Of course a 
particularly interesting case is mo = 0, which corresponds to a vanishing vector meson 
bare mass in the restored phase ( "Brown- Rho scenario"). Therefore, in the following we 
will present fits with m = as well as with the maximum possible value mo = 400 MeV. 
Together with the two values for the cr-meson mass this leads to four different parameter 
sets, which are listed in Table |I|. As explained in Sect. |3~T| , the scale parameter A can be 



chosen arbitrarily, since the dependence of the final results on A can be absorbed in the 
finite parts of the counterterms. On the other hand, the numerical values of the finite 
parts of the counterterms are meaningful only if the value of A is known, and this is the 
reason why it is listed in Table |I[ 
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Fit A 


Fit B 


Fit C 


Fit D 


Input 




m n / MeV 


140 


140 


140 


140 






m a 1 MeV 


600 


600 


800 


800 






U I MeV 


93 


93 


93 


93 








0.1 


0.1 


0.1 


0.1 






m / MeV 





400 





400 






A /MeV 


1000 


1000 


1000 


1000 


Obtained from 


m p / MeV 


767 


767 


767 


767 


fit of v\ and a\ 


m a , 1 MeV 


1064 


1062 


1056 


1054 


and from 




J pTTTT 


7.63 


7.77 


7.92 


7.97 


conditions 




(f+Jf) 


0.076 


0.078 


0.083 


0.084 


specified in the 




-9.36 


-9.30 


-8.48 


-8.56 


text of Sect. 


5.1 


6jj? / MeV 2 


6666 2 


3873 2 


6211 2 


3522 2 








-3147 


-1421 


-2671 


-1221 






JZy 


0.260 


0.064 


0.069 


-0.056 






{5ml + Jh 2 <?l) 1 MeV 2 


2352 2 


2197 2 


2245 2 


2128 2 






Sh 


20.1 


14.9 


15.4 


9.11 






Tg 


33.7 


32.5 


29.3 


29.0 


Computed from 


a / MeV 


150.9 


152.5 


153.7 


154.1 


parameters 




m w0 / MeV 


86.1 


85.3 


84.6 


84.4 


listed above 




A 2 


7.74 


7.58 


13.4 


13.3 






c/ MeV 3 


103.8 3 


103.5 3 


103.2 


103.1 






9 


5.55 


5.52 


5.47 


5.45 






h 


23.9 


23.2 


22.3 


22.0 








25.8 


18.4 


24.9 


18.0 






-iT/ml /MeV 


9.80 


14.40 


7.65 


10.82 



Table 1: Model parameters and related quantities in four different fits. 
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Figure 15: Left panel: the Real- and imaginary part of the p- meson propagator as a function 
of M = \fq* for the parameter sets A and D listed in Table El. Right panel: the tttt scattering 
phase shifts in the J = 1,1 = 1 channel for parameter set A (the three other parameter sets 



give the same result). The data are taken from Ref. [31] 



5.2 The p meson and the vector correlator 

The real- and imaginary parts of the p-meson propagators G p (q 2 ) = l/(g 2 — m 2 — X* (q 2 )) 
obtained for parameter sets A and D of Table [3] are displayed in the left panel of Fig. |l|. 
As a result of the fitting procedure described in the previous section the propagators look 
very similar for both sets. In fact, the curves for the two other parameter sets, B and C, 
lie in between. 

In the right panel of Fig. we show the tttt scattering phase shifts in the J — 1, 1 — 1 
channel. Near the resonance, s ~ m 2 , the tttt T-matrix in this channel is clearly dominated 
by the s-channel p-meson exchange. Thus the scattering phase shifts are given by tan<5j = 
ImGp/KeGp. As can be seen in the figure, the phase shifts obtained in this way fit the 
data up to 1 GeV. In this region the four parameter sets of Table p] lead to practically 
identical results. 

The electromagnetic form factor of the pion is not simply proportional to the p-meson 
propagator as in Sakurai's vector dominance model, in which the photon can couple to 
the two pions only via the p meson. Instead, we have to consider also the direct coupling 
of the photon to two pions. Thus the pion electromagnetic form factor is given by 



FM 2 ) = 1 + f^^T^ • (86) 



The result can be seen in the left panel of Fig. [H], together with the experimental data. 
Except for the small structure at q 2 = 0.61 GeV 2 caused by p — u mixing, the data are 
well reproduced up to q 2 ~ 1.5 GeV 2 (fit D) or even higher (fit A). 



In the right panel of Fig. [16] we display our results for the vector-isovector spectral 
function V\(q 2 ). In the resonance region, the experimental data are perfectly reproduced, 
but in the continuum region, above q 2 ~ 1 GeV 2 , we get only a small fraction of the 
strength measured by the ALEPH collaboration. This is obviously related to the fact 
that above 1 GeV four-pion channels become important. In our model, these are partially 
contained as irai, ap and pp channels. The threshold for the ap channel at q 2 = 1.87 GeV 2 
for fit A and q 2 = 2.46 GeV 2 for fit D is clearly visible, and also the pp channel gives 
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Figure 16: Left panel: the electromagnetic form factor of the pion for the parameter sets A and 
D listed in Table |. The curves for the parameter sets B and C lie in between. The data are 
taken from Refs. p2|, |33|. Right panel: the spectral function vi(q 2 ), together with data obtained 



from r decay by the ALEPH collaboration [22] 



a sizable contribution above 2.35 GeV 2 , whereas the nai channel contributes almost 
nothing. At first sight this is surprising, since in the p-meson spectral function the nai 
channel is relatively strong (in the left panel of Fig. [15] one can see the threshold as a 
small kink at 1.2 GeV). However, in the vector spectral function two amplitudes interfere 
destructively, namely (7 — > irai) and (7 — > p — > 71 a\). This interference is absent in 
the ap and pp channels, because the ap channel can be reached only via the p meson 
(7 — ► p — > ap), while for the pp channel only the direct coupling exists (7 — ► pp). 



5.3 The ai meson and the axial- vector correlator 



Before the precise r decay measurements by the ALEPH collaboration |22| , the knowledge 
about the axial vector correlator and the properties of the a\ meson was more or less 
uncertain, and even in the latest particle data book mass and width of the a% meson 



still have large uncertainties: m ai = 1230 ±40 MeV and T ai = 250 to 600 MeV. However, 
as pointed out in Ref. | 3"0[| , the discrepancy of different experiments to some extent stems 



from the use of different parametrizations of the a\ propagator in the analyses. Therefore, 
one must be careful when taking over parameters from another model. We have decided 
to fit the axial-vector spectral function ourselves and not to restrict the parameters to the 
values listed in the particle data book. 

The resulting a\ propagator G ai (q 2 ) = l/(q 2 — m 2 ai — E^g 2 )) is shown in the left 
panel of Fig. [17]. The mass, defined as the position of the maximum of — Im G ai , is rather 
low (see Table [I]). However, if we adopt another usual convention, defining the mass as 
the position of the zero of ReG ai , we obtain higher values between m ai = 1136 MeV (fit 
D) and m ai = 1159 MeV (fit A). The a\ meson width is defined as the difference of the 
two masses where -ImG ai crosses one half of its maximum. We find T ai = 424 MeV for 
fit A and r oi = 393 MeV for fit D. 

The spectral function of the transverse axial- vector correlator, ai(q 2 ), is displayed in 
the right panel of Fig. O. Between 0.9 and 2 GeV 2 we can fit the ALEPH data very well. 



In this range the spectral function is almost saturated by three-pion states f22| , which in 
our model are approximated by up states with sharp p mesons. Over a wide range of q 2 
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Figure 17: Left panel: the real- and imaginary part of the a\ propagator as a function of 
M = sfq* for the parameter sets A and D listed in Table [fl. The curves for the parameter sets B 
and C lie in between. Right panel: the transverse axial-vector spectral function a\ as a function 



of q . The data are taken from Ref. [22]. 



this seems to be adequate. Only in the region q 2 < (m n + m p ) 2 , of course, some strength is 
missed. Above 2 GeV 2 the spectral function receives sizeable contributions from five-pion 



states p2 |. These are partially contained in our model as aa± and pa\ states (for fit A, 
the aa\ threshold is visible at 2.77 GeV 2 ), but this is obviously not sufficient in order to 
fit the data. 

In Fig. [TS| it is shown that up to the r mass the exclusive three-pion final state can 
very well be approximated by a up final state with a sharp p meson. In the left panel we 
display the three-pion invariant mass distribution measured by the ARGUS collaboration 
|2l| and obtained in our model as the sum of 1/2 of the up and 2/3 of the 7ra spectrum 
(see App. 0). The na channel contributes only between 1.3 % (fit C) and 5.4% (fit B) to 
the integrated spectrum, which lies below the upper limit of 6 % given in Ref. p 



In the right panel of Fig. 18 we display the n n and 7r + 7r invariant mass distributions 
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Figure 18: Left panel: the three-pion invariant mass distribution of r~ — > tt^tt^it + ^ t for the 
parameter sets A and D listed in Table |l|. Right panel: the invariant mass distributions of tt^tt^ 
(theoretical curves and crosses) and 7r + 7r~ pairs (dots) of r~ — > ir~ir~iT + b> T . The data are taken 



from Refs. [21] and M. 
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Figure 19: The multi-particle contribution of the longitudinal axial-vector spectral function, 
a-o(Q 2 )i calculated with the parameters from fit A in Table ED. For the meaning of the different 
terms see Eq. (81). 



of the decay t~ — > tt~ty~ty + u t . The experimental data show that the ir + 7r~ distribution 
(dots) is the sum of the 7t~tt~ distribution (crosses) and a distribution which is concen- 
trated around the p-meson mass, indicating that the 7r~7r~7r + final state is indeed reached 
via a 7i~p° intermediate state. The interference term, which is responsible for the fact that 
at low invariant masses the 7r + 7r~ distribution lies below the 7t~tt~ distribution, is small. 
The 7r~7r~ distribution which corresponds to our model (see App. [B]) is in good agreement 
with the data. Since we do not take into account the p-meson width and the interference 
term, the corresponding 7r + 7r~ distribution would be given by the tt~tt~ distribution plus 
a delta function at the p-meson mass (and a small delta function at the cx-meson mass). 



5.4 The longitudinal axial-vector correlator 

As one can see from Eqs. (^) to (|84"D the longitudinal axial-vector correlator has a pole 
at q 2 = m^ 2 (= m 2 ) with residue Air 2 /^ 2 (= 47r 2 / 2 ), i.e. 

a (q 2 ) = Ait 2 f 2 J{q 2 - ml) + ~a (q 2 ) , (87) 



where ao(<? 2 ) denotes the contribution of multi-particle channels. According to Eq. (|8~TD, 
the longitudinal IF-boson self energy consists of three terms. The corresponding contri- 
butions to a (q 2 ) are shown in Fig. [T!| The first term gives a positive contribution to 
ao{q 2 ), whereas the other terms contain interferences of several amplitudes and are nega- 
tive. The sum of the three terms is of course positive, but almost zero. In fact, it is not 
visible in the figure since it is about a factor 10~ 4 smaller than the first contribution. To 
be more quantitative, we integrate the strength of a (q 2 ) up to m T = 1777 MeV. With 

2 

the parameters of fit A we find f™ T dq 2 do(q 2 ) = 75 MeV 2 , which has to be compared with 
the strength 4vr 2 / 2 = 0.34 GeV 2 of the single-pion contribution. 

The smallness of do(q 2 ) is a manifestation of the PCAC hypothesis. Usually the PCAC 
hypothesis is formulated as d^j^ = f^ml-n . However, beyond tree level this is almost 
meaningless, since the field operator ir a has non- vanishing matrix elements (/|7r a |0) not 
only for single-pion states |/), but also for multi-particle states, e.g. in our model for na, 
np, a\o and a±p states. In fact, as pointed out by Weinberg Ifffil , the PCAC assumption 
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should be interpreted in a different way, namely that do(q 2 ), which is directly related to 
the multi-particle content of c^j^jO), is negligible. Obviously this is fulfilled very well in 
our model. 

5.5 The (T-meson propagator 

A general problem of models with linearly realized chiral symmetry is that at tree level 
it predicts the existence of a sharp scalar isoscalar meson, which has not been found ex- 
perimentally. However, in the latest particle data book [^4| a meson with these quantum 



numbers is listed as "/o(400 — 1200) or cr". This particle has a very large width (ap- 
proximately as large as the mass) due to the decay into two pions. Since our calculation 
contains the a — > txix graph, one could expect that it should be possible to get a satis- 
factory description of the a meson. Unfortunately, this is not the case for the following 
reasons: 

The a7rir vertex has four contributions, analogous to the four contributions to the 
pirn vertex of the gauged linear a model shown in Fig. [I]. These amplitudes interfere 
destructively. If the outgoing pions are on the mass shell, the total vertex function has 
a zero at some value of the invariant mass M of the cr meson. For the parameter sets 
A and B listed in Table [l|, this zero lies at 754 MeV and 705 MeV, respectively. This is 
above the tree-level cr mass, m a = 600 MeV, but the amplitude is already very small for 
M pa m a . The authors of Ref. |19|] hoped that this effect could improve the description 
of the cr meson, since they believed that the width obtained from the linear cr model 
without vector mesons was too large. However, now the width is much too small, as can 



be seen in the left panel of Fig. |2(J, where we display the cr-meson propagator for the 
parameter sets A and B. The different normalization results from the fact that we adjust 
the wave-function renormalization constant 8Zq> for the pion propagator and not for the 
cr-meson propagator. 

Another problem arises from the fact that, at the two-pion threshold, the a- meson self 
energy usually becomes strongly attractive (i.e. negative). Since the cnnr vertex function 
at low q 2 is dominated by the ordinary a model vertex oc A 2 cr , with A 2 oc m 2 — m 2 , this 
effect is stronger for higher cr-meson masses and has therefore more dramatic consequences 
for the parameter sets C and D. For parameter set C we find a two-pion bound state at 
M = 231 MeV, while for parameter set D the cr-meson propagator has the wrong sign 
at low M and a pole in the complex plane. For illustration, we have plotted the inverse 
cr-meson propagator in the right panel of Fig. |2(J 

Obviously the cr-meson cannot be described in one-loop approximation. A good de- 
scription of the cr-meson in the linear a model can be obtained by summing the 7T7t loops in 
the self energy to all orders (RPA) in the large- N limit (N = number of pions) ||36|| . Sum- 
ming the 7T7T loops reduces the strong attraction in the real part at the two-pion threshold. 
However, the application of the RPA to our Lagrangian would be extremely difficult. On 
the other hand, for the vector mesons, in which we are interested, the one-loop calculation 
seems to be sufficient. 



6 Conclusions 

We have presented a chirally symmetric hadronic model for the vector and axial-vector 
correlators in vacuum. The starting point was the linear cr model extended by p-meson 
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Figure 20: Left panel: the real- and imaginary part of the cr-meson propagator for parameter 
sets A and B. Right panel: the real- and imaginary parts of the inverse cr-meson propagator for 
parameter sets C and D. 



and cti fields in a chirally symmetric way. Unlike the gauged linear a model, which leads 
to phenomeno logical difficulties (see Sect. our Lagrangian is only symmetric under 
global, but not under local SU(2)l x SU(2)r transformations. After all, chiral symmetry 
is also only an (approximate) global symmetry in QCD. 

For a realistic description of p and a± we could not stay at tree-level, where all mesons 
are stable particles. In order to include the most important decay channels, p — > tttt and 
a\ — > pn, we performed a one-loop calculation. Here care had to be exercised, not to 
destroy the symmetries of the model by choosing an unsuited regularization scheme. We 
have decided to use dimensional regularization in order to isolate the divergences, which 
are then subtracted by chirally symmetric counterterms. 

The price to pay for giving up the requirement of a local chiral symmetry is the loss 
of the so-called current-field identity which directly relates the vector and axial-vector 
currents to the meson fields. Therefore some additional work was needed to construct the 
correlators in a consistent way. 

In the last part of the paper we have presented numerical results. We have fixed our 
parameters to reproduce the empirical values of m n and f n and the peak regions of the 
vector and axial vector spectral functions, which have been extracted by the ALEPH col- 
laboration from the r decay. Basically, this leaves two parameters undetermined, namely 
the mass of the cr-meson, which is not a well established resonance in nature, and the 
mass parameter m . The latter is the fraction of the bare p-meson mass which is not 
proportional to the vacuum expectation value of the a field, i.e. the fraction which would 
survive in the restored phase. From our fit we only found an upper limit of 400 MeV 
for m . On the other hand, a particularly interesting case is mo = 0, which corresponds 
to a vanishing bare vector meson mass in the restored phase ( "Brown- Rho scenario"). 
Therefore, we have performed calculations with ttlq = and with mo = 400 MeV, both 
for two different cr-meson masses, m a = 600 MeV and m a = 800 MeV. 

All fits describe the peak regions of the experimental spectral functions equally well 
and also give a satisfactory description of the pion electromagnetic form factor, the tttt 
phase shifts in the vector-isovector channel, and the 7r~7r~7r + and 7r~7r~ mass distributions 
for the decay r~ — > ix~ n~ ix + v T . At higher energies, q 2 > 1 GeV 2 in the vector channel and 
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q 2 > 2 GeV 2 in the transverse axial-vector channel our fits underestimate the empirical 
spectral functions. In these regions four-pion (vector channel) or five-pion (axial-vector 
channel) final states become relevant, which are not properly contained in our one-loop 
calculation. 

The longitudinal axial correlator is completely dominated by the pion pole, in agree- 
ment with the PCAC hypothesis. The remaining strength, which is due to the small 
explicit breaking of chiral symmetry in our Lagrangian, is the result of a delicate can- 
cellation of three terms which are four orders of magnitude larger than their sum. This 
nicely demonstrates the consistency of our model with chiral symmetry as well as the 
stability of the numerics. 

In the present paper we have restricted ourselves to vacuum properties. However, 
having the technology developed, it is now straightforward to calculate the correlators at 
finite temperatures. Work in this direction is in progress. The extension to finite baryon 
density in a chirally symmetric way is obviously more difficult. 
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A The three-pion invariant mass distribution of r — ► 

The differential r decay width has the form 

= G ^ 1 ife^("^-"+ + « + ™i>*— ) ■ m 

In this formula, Gf = e 2 / (4\/2 sin 2 Owm-w) denotes Fermi's coupling constant and H l n - n _ n+ 
and if£ + are the longitudinal and transverse components of the so-called hadron ten- 
sor, respectively. Under the assumption of PCAC, the longitudinal component can be 
neglected. 

In our model the 7r~7r~7r + spectrum consists of 1/2 of the np plus 2/3 of the na 
spectrum, i.e. 

Ht m+ {Ml) = l -Hl p {Ml) + \hUmI) . (89) 
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The up and ira hadron tensors are given by 



KM 1 ) 



1 f QnQti' _ 

' 3 ^ n 2 w 



d 1 f d kp 1 



r - /j(2 7 r)3 2A;07(2^2io (27r)4 ^-^-^ 

16c 



X Z„i[ PV P!/ - ) (y4 l4 / 7rp )^ c (A^ p ) 



^ 2 



1 / g M g^ 
3 V - 2 



g - " J J (2vr) 3 2fc°y (2tt) 3 2&° 



(27r) 4 <5(g - K - K) 



with ki 



sin 6w 
e cos 6c 
sin # 



m 2 + k 2 and 

( r ^p)a6c + ( r ^i)a^( Ga i(^)) K A( ra i^) 



A;/ 
dbc 

db 



The vertex functions Twn P and r^Tro- are defined analogously to Fig. 0. 
The curves shown in the left panel of Fig. 113 are normalized as follows: 



dN„ 



ir~ir~'n-< 



dM- 



N T -T T - 



dr 



N, 



dr 



37T 



3tt 



B( 



71 7T 7T^^ 



3tt 



(90) 



(91) 



(92) 
(93) 



(94) 



7T 1T 7T T 



37170 is the total number of n it h + events ||21||, r T - = 290.0 fs is the 



where N. 

t lifetime [2^|, and _B(r~ — > 7r~7r~7r + z/ T ) = 9.15% is the branching ratio |22[ . 



B The 7r 7r invariant mass distribution of r — > 

7T~7T~7T + h f T 

Since our model does not contain the 7r~7r~7T + final state, but only 7T~p° and i\~o final 
states, the t\~ix~ distribution cannot be computed directly. Instead we must interprete 
the p° or a meson in the final state, with momentum k p or k a , respectively, as a 7r + 7r~ 
pair with momenta p n + and p w - (i.e. p n + + p w - = k p or k a , respectively). For simplicity 
let us first look at the ira final state. 

In the rest frame of the a meson the kinematics is given by 



k T: 



(m a , 0) , 

(q -m a ,q), q 



i i i 
q — m a — m^ 



a/ (q 2 — m\ — ml) 2 — ^.m\m 2 a 



2m r . 



Thus the 7i~7i~ invariant mass depends only on q 2 and z = cos(]9, q): 

2 q 2 - ml + 3ml 



(95) 



M. 



2tt 



(96) 
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Let P a {z) be the probability for a certain angle between p and q. Since the a meson 
decays into a 7t + 7t~ pair in s wave, the angular distribution in its rest frame is isotropic, 
i.e. P a (z) = const. = 1/2. Now the tt<j contribution to the 7t~tt~ invariant mass spectrum 
can be written as 

dN„ f 00 , 2 /dNj:Z>, 



dMl J 9ml V dMl J vl=i 2 

dzP a {z)S(Ml - (q 2 -ml + 3ml) /2 + 2\p\\q\z) 



(97) 

where dN^.^/dM^ is the ttu contribution to the three-pion spectrum discussed in 
App. [A], and the integration boundaries q± in the last line are given by 

2 2m^ + mjMl ± 2\p\m a M2n\/Ml - Ami 
«± = h^2 • (98) 

For the np contribution the situation is very similar if one replaces everywhere m a by 
m p . The only complication comes from the fact that the angular distribution cannot be 
assumed to be isotropic, since the p meson decays into a 7r + 7r~ pair in p wave. In order to 
obtain the angular distribution, we assume that the p meson is measured by a fictitious 
detector, which forces it onto its mass shell without changing its spin, and the decay 
into two pions takes place after this measurement. Under these conditions the angular 
distribution is given by the square of the amplitude 

B^e abc = (A Wwp )^(p w - „ - p w + „) , (99) 

averaged over the three transverse W polarizations, and evaluated for the kinematics 
given above (but with m a replaced by m p ): 

P p {z) oc (M£ - gJjBW* = c + c 2 z 2 . (100) 



The normalized angular distribution is then given by 

c + c 2 z 2 
2(c + c 2 /3) ' 



P P (z) = J° + C2Z l, • (101) 
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